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ABSTRACT 



Numerical wave refraction programs permit a detailed study of the 
transformation of wave energy as waves move from deep water to shallow 
water. By eliminating the subjectivity that is present with hand drawn 
diagrams the effect of small variations in the initial assumptions and 
wave conditions can be investigated. The effects on wave refraction of 
tide level changes and bottom friction are investigated here. It is 
demonstrated that a uniform increase in water level as would be caused 
by tidal fluctuations can cause a significant change in the wave re- 
fraction pattern for a given nearshore region. A computing procedure 
is developed to permit numerical refraction programs to account for 
bottom friction. ■ The reduction in wave height caused by bottom friction 
depends primarily on bottom slope; this relation is shown in tabular 



and graphical form. 
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Introduction 



Preparing accurate wave refraction diagrams is a necessary first 
step in many coastal engineering and amphibious operations. The advent 
of the numerical wave refraction programs [[Griswold (1963) , Stouppe 
(1966), Wilson (1966)J has considerably reduced the time required to 
construct a wave refraction diagram. The speed of the computer opera- 
tions permits more parameters to be included in the calculations and 
improves the accuracy of the output. Of equal importance is the fact 
that numerical products provide a completely objective basis for evalua- 
ting the cumulative effect of various assumptions on wave refraction. 

The present study has investigated the effects of tide-level 
changes and bottom friction on wave -re fraction patterns for several 
environments. The investigation first compared the wave refraction 
programs of Griswold (1963) , Stouppe (1966) , and Wilson (1966) . Since 
there was no way available to the author to compare program accuracy 
on an absolute scale, the comparison was based on method of computation, 
simplicity of use, and presentation of results. 

The program written by Stouppe was considered to have the most 
potential and was chosen for use in this investigation. The major 
advantages of Stouppe *s program are these: 

(1) The program uses water depth directly as the interpolation 
surface for evaluation of wave speed, thereby limiting the 
number of computations and increasing the program versatility. 

(2) The program uses a second-order non-linear differential equation 
to determine the refraction coefficient at each point along 

the ray. Thus, the wave height at each point can be computed. 
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(3) The program uses constant time steps in the computation rather 

j than constant distance steps. This reduces the distance be- 

tween computations as the shoreline is approached. Thus for 
the same number of computer operations a greater percentage of 
computations are in the shallowest area where refraction, 
friction, etc., are most important, and are changing rapidly. 

(4) The program includes a graphical output that plots the wave 
orthogonals, wavecrests and the shoreline. 

The following modifications were made to Stouppe's program: 

(1) The program was made to recycle, permitting more effective use 

of computer time when more than one diagram is to be constructed. 

(2) Provision was made to add a constant value to the depth field 
to investigate tide level effects. 

(3) A subroutine was added to account for bottom friction in the 
computation of wave height. 
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2, Computer Program, 

Depth grid . 

The first step in the utilization of this program is the construction 
of a grid of water-depth values for the desired area. The grid must be 
sufficiently large so that the starting point of all rays to be followed 
is in deep water, i.e., the ratio of water depth to deep water wave 
length (d/Lo) is greater than 0,5, By convention the x-axis is positive 
increasing toward the shore while the y-axis is positive to the left of 
the x-axis. The grid interval is selected so that the bottom contours 
are reasonably parallel to one another within a given square. 

In constructing the depth grid it is advantageous to record the 
average value of all depths within that grid square, rather than the 
value of water depth at the intersection point. This procedure tends to 
correct for small variations in the depth of water within the depth grid 
without necessitating the use of a smaller grid size. All actual depths 
are made positive; extrapolated depth values are continued on land for 
two grid units from the shoreline and are made negative. Beyond two 
grid units from the shoreline any arbitrary negative depth may be used. 

For depths on the shoreline itself zero is used. 

To reduce the time required to obtain the depth grid a clear plastic 
overlay with a black dot at each grid intersection point was prepared by 
initially locating the grid dots on white paper and producing a viewfoil 
transparency with an Ozalid copier, Model 400. 

Input data , 

Stouppe's order of data input was modified to permit the program 
to recycle. The present data deck consists of three sections. The first 
section contains only the value of MM and NN which designate the number 
of points in the depth grid in the x and y directions respectively. The 
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second section contains the water depth grid. These values are arranged 
so that the computer reads all the y values for each successive x-grid 
position. Water depths are read in fathoms and are converted to feet 



within the computer program. The last section contains one card for 
each set of input parameters for which it is required to construct a 
wave refraction diagram. The input parameters include: 



X,Y: 


The grid position for the starting point of the first ray 
(feet) . All starting rays must be more than two grid 
intervals in distance from the edge of the depth grid. 


FCF: 


The coefficient of bottom friction. 


TIDE: 


The tide level, in feet, above chart depth for which com- 
putation is required. This value is subtracted from the 
depth grid after each refraction pattern is completed. 


! HINT: 


The deep water wave height (feet) . This is required for 
the friction subroutine. 


T: 


The initial wave period (seconds) . 


Al: 


The initial wave angle (degrees) . The smaller angle 
measured between the positive Y direction and the initial 
wave crest. The angle is made negative if the slope of 
the wave crest is negative. 


NOR: 


The number of rays to be followed. 


TIME: 


The time interval between computations for advancing of 
the wave front (seconds) . 


DIST: 


The distance between rays (feet) . 


GRID: 


The grid interval (feet) . 



This order of data input permits wave refraction diagrams to be 
computed for many combinations of input variables while requiring the 
computer- to assemble the program and read the depth field only once. 
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G omput e r out put , 

The output from the computer consists of both a printed output and 
a graphical output. The printed output variables are these: 

X,Y: The x and y coordinates for a given wave crest and ray 

number (yards) . 

COREFR: The coefficient of refraction (K^) for the wave at the 

point x,y. 

HHO: The ratio of wave height to deep water wave height 

neglecting bottom friction. 

HHFR: The ratio of wave height to deep water wave height when 

bottom friction is included. 

NGO: Indicates that the ray has terminated (1) , or it is 

continuing (2) . 

DEPTH: The water depth for grid position X and Y (feet) . 

The graphical plot of wave crests is programmed for the Calcomp 
160 system utilizing the DRAW subroutine programmed by J. R. Ward for 
the Fortran 60 system. See Appendix II, Stouppe (1966) , for subroutine 
listing. The first and each succeeding third wave crest are plotted. 
The zero water depth points are used for contouring the shoreline on 
the graph. 

Appendix 2 contains a sample of the printed and graphical output 
for Monterey Bay, California. 

Computer operations . 

Since the details of the computer operations are discussed by 
Stouppe (1966), only the main features will be discussed here. The 
formulas used in the computations are listed in Appendix 3. 
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The water depth at the first point is computed by fitting the 
closest nine grid-point depths to a quadric surface by the least-squares 
method. An iterative procedure is then used to solve for wave velocity. 
It is assumed here that wave velocity is a function of water depth and 
period only. 

The wave ray is then moved to the next point by solving for the ray 
curvature and projecting the ray forward a distance equal to the product 
of the wave speed for the point multiplied by TIME, At this new point 
the value of ray separation (^), coefficient of refraction (K^) , 
shoaling factor (Kg), H/Hq, and HHFR are calculated. 

This procedure is repeated until all orthogonals have been advanced 
one time interval. The new wave crest is plotted and the entire pro- 
cedure is repeated until all orthogonals have either gone off the grid 
or reached the shoreline. 
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3. The Effects of Tide Level Changes on Wave Refraction Patterns. 
Differences in shoaling factor , wave speedy and ray curvature at specific 
depths . 

The effect on wave refraction of a uniform increase in water depth 
caused by tide level changes is cumulative and is dependent on the bottom 
depth contours. Therefore, to evaluate the effect properly, a wave ray 
must be followed along its entire distance from deep water to the shore. 
However, by calculating the value of shoaling factor (Kg), wave speed (C) , 
and ray curvature (FK) at both high tide and low tide it is possible to 
evaluate the order of magnitude of the differences in the basic parameters 
of wave refraction caused by tide level fluctuations. 

The value of Kg, C, and FK were calculated initially for 25, 50, 
and 100 foot water depths; second calculations were made for a water 
depth equal to the initial water depth plus a tide level increase of 
6, 10, and 16 feet. To evaluate these differences on a relative scale 
a percent difference was calculated by using the formula: 



% difference X = 



X H r ^ LT 

X H T 



( 1 ) 



where 

X = variable representing C, Kg , or FK 
Xj^t = value of X at high tide water depth 
Xxt = value of X at low tide water depth. 

The percent difference of Kg was generally small and under certain 
conditions went to zero (for those values associated with the inflection 
point in the curve when Kg is plotted against water depth) . 

Figure 1 shows the percent difference in wave speed corresponding 
to the low tide water depth for a specific tide range and wave period. 
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Figure 1. EFFECT ON WAVE SPEED OF VARYING TIDE LEVELS FOR 8-16 SECOND WAVES. 
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WATER DEPTH (LOW TIDE) 



It is apparent that the difference is often greater than 10 percent and 
increases with increasing wave period, increasing tide range and decreasing 
low tide water depth. 

The values of ray curvature (FK) were determined using the equation: 



F K - i C'^F) 



AC. _ 



/oc<jQ dA. 



7 



( 2 ) 



where 



A = approach angle 

4^ . = derivative of wave speed in the X and Y directions 

respectively. 

Figure 2 shows the percent difference in FK versus tide range (computed 

^ c 

using equation (2) with the following initial values: A = 30°, = 0, 

= A C for a twenty foot change in water depth starting with an 
initial depth of 100 feet) . The values of percent difference plotted 
are representative of the percent difference for all values of 
A, j ^ j and water depth: this occurs because the terms inside 

the brackets in equation (2) have nearly the same value for high tide 
and low tide and tend to cancel out in the percent difference calculation. 
This reduces the percent difference calculation for FK to 



% difference FK = 




/■* 



(3) 



^ nr 

The fact that the difference calculation is primarily dependent on tide 
range and wave speed is clearly shown in Figure 2. 

It is concluded that the change between high tide and low tide in a 
refraction pattern is primarily determined by the change in ray curva- 
ture (FK) , 
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Figure 2. EFFECT ON RAY CURVATURE OF VARYING TIDE LEVELS FOR 8-16 SECOND WAVES 
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TIDE RANGE (FEET) 



Difference in wave refraction patterns for specific locations . 



To investigate the combined effect of the changes in C, Kg and FK 

caused by a uniform increase in water depth, wave refraction diagrams 

were constructed for several locations using both the low tide water 

depths and the high tide water depths. 

The first area considered was the California coast near Big Sur, 

between latitude 36° 14* 30” and latitude 36° 17* 42”. This area was 

chosen because eight -second waves arriving on the coast will be refracted 

by the offshore bar (shown by the shaded area in Figures 3 and 4) at low 

tide but not at high tide.. Thus tide level changes will cause selective 

refraction by the offshore bar. Figure 3 shows the refraction pattern 

for 8 second waves approaching from 270 degrees at low tide; Figure 4 

shows a different pattern for the same condition at high tide.^ 

Figures 3 and 4 show two major differences in the refraction pattern. 

First, orthogonals 6 and 7 cross approximately halfway to shore at low 

2 

tide but they do not cross at high tide. Secondly, at low tide there 
is a convergence of wave orthogonals at Point A; this convergence is 
not present at high tide. 



^To make the changes in the orthogonal pattern more prominent 
Figures 3-8 are plotted with an expanded ordinate scale. Also, the 
computer was not used to contour the wave crests but rather the wave 
orthogonals were contoured manually; the shoreline was indicated by a 
triangle plotted whenever water depth was zero. TIME equalled 30 
seconds for all diagrams. 

2 

Interpretation of the effect of crossed orthogonals is beyond the 
scope of this paper. In this discussion we will only be concerned with 
the difference in ray paths between low tide and high tide. 
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Table 1 shows the value of the wave-height ratio (H/Hq) computed 
at four wave crests before the shoreline is reached for specific points 
along the coast. This table clearly shows that the value of wave height 
can be considerably different at low tide and high tide. 

Figures 5 and 6 represent the same initial conditions as Figures 
3 and 4 except that the initial wave direction is from 310 degrees. 

The change in the orthogonal pattern is clear. It is noted that in 
Figures 3 and 4 the primary change in the orthogonal pattern (between 
low tide and high tide) is in the vicinity of points B and C, whereas 
in Figures 5 and 6 the greater change occurs near Point A. 

Wave refraction diagrams were also computed for the coast of Maine 
between latitude 44® 39^ 35'* and latitude 44® 44^ 45*'. This area was 
chosen in order to determine if the combination of a large tide range 
(16 feet) and deep water depths (greater than 200 feet a short distance 
from the shoreline) would cause a change in the refraction pattern. 
Figures 7 and 8 show the refraction patterns for 16-second waves arriving 
from 070 degrees at low tide and high tide. The difference in refraction 
patterns is evident. Table 3 shows the H/Hq ratios (calculated at four 
wave crests before reaching the shoreline) for specific locations at 
low-tide and high-tide conditions. 

In total, about 20 sets of refraction patterns were analyzed and 
it appears impossible to predict what changes in wave height and re- 
fraction patterns will occur as the result of tide level changes. The 
only generalization that can be made is that the curvature of an or- 
thogonal is generally greater at low tide than at high tide. It is 
concluded that precise water-level information, including tide state, 
should be incorporated into wave refraction procedures. 
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TABLE I 



H/Hq ratio at low tide and high tide for specific coastal locations 

NEAR POINT SUR, CALIFORNIA. 



LOCATION 
(SEE FIG. 3) 

A 

B 



H/Ho RATIO 

LOW TIDE HIGH TIDE 



1.1 

1.3 



1.2 

.9 



C 3.1 1.7 

D .8 .8 



TABLE II 

H/Hq ratio at low tide and high tide for specific COASTAL LOCATIONS 
NEAR EASTERN HEAD, MAINE. 



LOCATION H/Hq ratio 



(SEE FIG. 7) 


LOW TIDE 


HIGH TIDE 


A 


2.8 


3.3 


B 


1.0 


CM 

• 


C 


,7 


.8 


D 


.7 


.7 



23 




x: 

o 

c 

•H 



• 4J 
U fJU 
<l) 

OT O 

o 
o o 

CO 

II O) 

»— I 
u CD 

§ o 
M W 
H 

x: 



o 


Z3 


o 


O 




C/!) 


CvJ 


1 






II 


4J 






< 


O 






• 




o 


• 


0) 


rC 


CO 


o 




c 


00 


•rH 


II 






4-J 


H 




. 


o 


CD 


o 


•rH 


o 


c 


vD 


u 




o 


Q) 


CM 


T— 1 


•rH 


CD 


t— 1 


O 


CD 


C/3 


o 






4J 


•N 


CO 




0) 




:2 


W 


1 




4J 


00 


CO 


•rH 


CD 


pq 


W 




• 


O 




M-l 


CvJ 




•<^ 


B 




CD 




J-4 


»— 1 


ooo 


CD 


vO 


•rH 


CO 


"O 






<U 


c 


"O 


o 


ZJ 


•rH 


4J 


4J 


•rH 


U 


4J 


CD 


CD 


U 


»— 1 


M-l 




0) 


O 




4-J 


0) 


- 


'O 


o 


•rH 


CO 


4-J 


— 




vl* 




»— 1 


o 


o 


hJ 


VO 




CO 


CO 


0) 




TJ 


<L) 


0 




4-J 




•rH 


00 


4-J 


•rH 


CD 







24 




23 



Figure4. High tide refraction diagram for Big Sur, California. T = 8 sec. A = 270°. TIME = 30 sec. 
Latitude 36°14'30" to 36°17'42". East-West scale 6000 Ft. /inch. North-South scale 4000 Ft. /inch. 
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Figure 5. Low tide refraction diagram for Big Sur, California. T = 8 sec. A = 310°. TIME = 30 sec 
Latitude 36°14'30" to 36°17'42". East-West scale 6000 Ft. /inch. North-South scale 4000 Ft. /inch. 
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Figure 6. High tide refraction diagram for Big Sur, California. T = 8 sec. A = 310 . TIME - 30 sec. 
Latitude 36°14'30” to latitude 36°17'42". East-West scale 6000 Ft. /inch. North-South scale 4000 Ft. /inch. 
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Figure?. Low tide refraction diagram for Eastern Head, Maine. T = 16 sec. A = 070°. TIME = 30 sec. 
Latitude 44°39'35" to latitude 44°44'45". East-West scale 8000 Ft. /inch. North-South scale 6000 Ft. /inch. 




29 



Figure 8. High tide refraction diagram for Eastern Head, Maine. T = 16 sec. A = 070^. TIME = 30 sec. 
Latitude 44^39*35” to latitude 44^44*45**. East-West scale 8000 Ft. /inch. North-South scale 6000 Ft. /inch. 



4. Bottom Friction Effects on Wave Height. 

Introduction to bottom friction . 

The present method for calculating wave height assumes that the 
propagation of energy between two adjacent orthogonals is conserved as 
the wave proceeds from deep water to shallow water. This implies that 
the addition of energy by surface winds, dissipation of energy by bottom 
friction or internal friction, and the reflection of energy by a steep 
slope are absent. In nature this is not the case - the wave energy 
between adjacent orthogonals may be changed by all of these means. 

According to Lamb (1932) , in fluids of small viscosity (such as 
sea water) , the velocity of the wave and the relationship between depth, 
length, and period are unaltered by internal friction. This is based on 
the assumption of a sinusoidal wave form but the result appears to in- 
dicate a very slow rate of damping of all waves in deep water, which is 
in accord with the limited observations in nature (Putnam and Johnson, 
1949) . 

At the sea bottom, however, friction can cause a loss of wave energy. 
It can be demonstrated both theoretically and experimentally that the 
particle velocities of shallow water waves are finite at the bottom and 
actually approach the theoretical orbital values at a short distance from 
the bottom. Therefore, a boundary layer in which mechanical energy is 
dissipated through the distortion of fluid elements must exist at the 
bottom. The effective roughness of a sandy bottom probably is determined 
by the size of the ripples in the sand rather than by the size of the sand 
grains. The ripple size, however, appears to be determined by the size 
and sorting coefficients of the sand in addition to the height and period 
of the wave and the depth of water. 
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Putnam and Johnson established the general magnitude for the dimen- 
sionless friction coefficient for the sea bottom to be 0.01. (This value 
was computed assuming initial conditions of 12 second waves, 4 feet high, 
traveling over a sandy beach with a ripple spacing of 5 inches.) Putnam 
and Johnson point out that since the friction factor only enters the 
equations for wave height to the first power, the friction factor 
could change by + 50 percent without appreciably changing the order of 
magnitude of wave-height reduction. 

Based on the following assumptions Putnam and Johnson were able to 
develop a procedure to calculate wave height when bottom friction is 
included . 

(1) The oscillatory motion at the bottom is sinusoidal and 
corresponds to the wave theory based on small amplitudes. 

(2) The friction coefficient does not vary with velocity and 

is approximately independent of water depth. This assumption 
amounts to choosing an average friction factor for the entire 
region over which the sea bottom affects the wave motion, 

(3) Perpendicular flow at the bottom resulting from percolation 
in the bottom material is negligibly small, 

(4) The sea bottom has a constant upward slope from offshore to 
the breaker point. 

Computing procedures . 

Bretschneider and Reid (1954) extended the original work of 
Putnam and Johnson and developed new equations to solve for wave 
height which included simplified equations for the special cases of 
constant water depth and constant bottom slope. 
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Bottom friction is of importance only on gently sloping beaches 
and in shallow water; thus it is possible, when using numerical wave 
refraction procedures, to calculate the reduction to wave height due 
to bottom friction by making the assumption that the water depth is 
constant between successive computations along the ray path from deep 
water to the shoreline. 

"'‘‘The wave height at each step is computed by the equation 



H = H, [ I- [d-KfWO- K4p-» ■ • • --0-K4J}] (4) 

where 



= the calculated wave height, neglecting bottom friction, for 
the point; 

Kf^ = the refraction factor between point n-1 and n. 

The coefficient of friction (K^) is calculated by the equation (from 
Bretschneider and Reid) 



Kf 








•f H. 

T" 



9 



( Ks 



(5) 



Kg - direct shoaling factor; 

= incremented distance. 

Kf is computed at each ray point with Ax equal to the distance 
along the wave ray between point n and n + 1. It is assumed in the 
computations that the water depth between point n and n + 1 is constant 
and equal to the depth at point n. 



32 



Results of bottom friction computations. 



The accuracy of the approximation that water depth is constant can 
be seen by a comparison of wave height obtained by this procedure to 
the wave height obtained by the more sophisticated and accurate pro- 
cedure used by Putnam and Johnson. This comparison is shown in Table 3 
(computed for wave period = 12 seconds, bottom slope = 5 feet). The 
comparison is very favorable along the entire distance of travel from 
the point that the wave first feels bottom to the breaking point. The 
difference is certainly within the order of magnitude of the error in- 
troduced by other assumptions. The difference in values is caused by 
the approximation to constant depth steps in the calculating procedure. 

Since this program uses constant time steps in moving each ray point 
toward the shoreline, the horizontal distance over which depth is assumed 
constant can be reduced by shortening the computing time interval (TIME). 
Table 4 shows the reduction (in percent) of wave height at the breaking 
point for different computing time intervals (with wave and bottom con- 
ditions the same as shown for Table 3 above). As expected, there is a 
slight increase in wave height reduction as the computing interval is 
decreased. It can be concluded that the value of TIME will affect the 
wave heights calculated by this procedure, but the effect is small. 

Realizing the limitations imposed by the initial assumptions of 
Putnam and Johnson and the calculating procedure employed here, the 
computer program was used to compute the wave height (with bottom friction 
included) for varying initial periods and bottom slopes. The results are 
shown in Table 5. (In arriving at these results it was assumed that the 
initial wave height was 5 feet and the friction coefficient was 0.01.) 

The reduction in wave height (in terms of percent of the wave height 
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TABLE III 



WAVE HEIGHT REDUCTION DUE TO BOTTOM FRICTION AS DETERMINED BY PUTNAM 
AND JOHNS'ON AND BY THE NUMERICAL WAVE REFRACTION PROGRAM. 



T = 12 sec. Hq 


= 5 FT. f = .01 


TIME 


= 30 sec. j 


Bottom slope = 1:300 


Breaking depth = 


11.5 


Ft. 


DISTANCE FROM 
WHERE h = Lo/2 
TO BREAKING POINT 


PERCENT REDUCTION 
FARRELL 

CONSTANT DEPTH STEPS 


WAVE HEIGHT 
PUTNAM & JOHNSON 
CONSTANT SLOPE 


PCT 


PCT 




PCT 


0.0 


0.0 




0.0 ' . 


50.0 


0.0 




0.0 


62.0 


0.2 




0.2 


70.3 


0.4 




0.4 


80.5 


1.0 




1.1 


86.8 


1.9 




2.1 


90.8 


3.1 




3.2 


93.0 


4.2 




4.5 


95.0 


6.0 




6.2 


97.1 


8.7 




9.3 


100.0 


18.5 




20.8 ' 



! I 



34 



TABLE IV 



PERCENT REDUCTION IN WAVE HEIGHT FOR VARIOUS VALUES OF TIME 
TIME (SEC) WAVE HT. REDUCTION ( %) 

20 19.0 
30 18.5 
40 18.0 
50 17.0 



TABLE V 

WAVE HEIGHT REDUCTION FOR BOTTOM SLOPES OF 1:100 to 1:400 AND WAVE 



PERIODS OF 8 to 


16 SECONDS 


(PERCENT) 






WAVE 




BOTTOM 


SLOPE 




PERIOD 


1:100 


1:200 


1:300 


1:400 


8 


5.4 


14.0 


21.5 


28.6 


10 


5.2 


13.5 


20.2 


27.0 


12 


5.2 


13.0 


18.5 


24.8 


16 


5.0 


12.2 


18.2 


24.2 


AVERAGE 


5.2 


13.2 


19.6 


26.9 
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without friction) is shown for the theoretical breaking depth (from 
H.O. 234). To obtain the percent reduction at the specific breaking 
depth, it was necessary to fit a smooth curve to the plotted points of 
wave height reduction at the three closest water depths and interpolate 
the wave height reduction at the breaking point. 

Figure 9 shows the average reduction (for all periods from Table 5) 
of wave height versus bottom slope. The curve is extended by a dashed 
line to the reduction amount found by Putnam and Johnson for a 1:10 
beach slope. 

The quantitative value of Figure 9 is controlled by the input 
accuracy of the coefficient of friction, the initial assumptions and 
the computation procedure. However, Figure 9 does show, qualitatively, 
the effect of bottom slope on wave height reduction due to bottom 
friction. Figure 9 is considered sufficiently accurate to be used in 
surf forecasting on sandy beaches to provide a first approximation of 
wave height reduction due to bottom friction. 
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PERCENT REDUCTION 



Figure 9 . EXPECTED WAVE HEIGHT REDUCTION AT BREAKING 
POINT DUE TO BOTTOM FRICTION FOR BOTTOM SLOPES OF 
1:10 to 1:400. 




^7 



5 . 



Conclusions . 



The numerical approach to wave refraction is a powerful tool in 
that it eliminates human subjectivity and permits the refraction com- 
putations to account for more factors than is possible with hand drawn 
techniques. 

A computer program that uses constant time intervals between com- 
putations has more suitability in accounting for the effects of friction, 
wind, reflection, etc., than a computer program that employs constant 
distance steps. 

Tide level changes can have a significant effect on wave refraction. 
This is especially true in areas where there is a large tide range or 
an offshore bar that can cause selective refraction. 

The changes shown by the tide level effects also point to the im- 
portance of having an accurate water depth grid. To avoid significant 
errors in the refraction pattern a great deal of care must be exercised 
in the preparation of the depth grid. 

Bottom friction can and does have an effect on the wave height. 

The numerical refraction program provides a convenient way to account 
for this reduction. 
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6. Recommendations. 

The course of this investigation has permitted a close look at wave 
refraction from many points of view. The numerical refraction program 
is capable of yielding wave refraction information with great speed and 
accuracy. (Computer time is approximately 5 minutes for a 20 ray diagram.) 
However, before the advantages of the computer product can be fully 
realized the author believes further investigations and improvements are 
needed in the areas of input data, computing procedures and field in- 
vestigations . 

The need for accurate depth information is essential to the compu- 
tation of accurate refraction patterns. It is anticipated that the present 
efforts of the Naval Oceanographic Office and the Coast and Geodetic Survey 
in the area of color photography will improve the water depth information 
especially in areas where the bottom contours change considerably and 
in the shallow areas where it is difficult for hydrographic vessels to 
proceed. 

The computer program can be expanded to account for more variables 
and to provide additional output information^ in particular: 

(1) The program should account for the fact that in nature there 
usually exists a spectrum of wave periods and directions, 
rather than the single period and direction that is considered 
here. 

(2) For amphibious operations it would be advantageous if the 
program determined the breaking depth of the waves, printed out 
the height of the waves at the breaking point, and then contoured 
the breaking point on the output refraction diagram. 
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(3) Another desirable output for amphibious operations is a computer 
plot on the output diagram of the area where the wave height is 
a minimum, as this would be an ideal area to anchor large ships. 

(4) ' The various subroutines should be checked to ensure that each 

section is providing comparable accuracy. 

Field work is necessary to determine the value of the bottom friction 
factor for various type bottoms. This is essential before accurate wave 
height information can be obtained. Field investigations are also needed 
to determine the accuracy of the computer program; until now the numerical 
wave refraction procedures have been checked only by comparison with hand 
drawn results. 



} 
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APPENDIX 1 



Computer Program for Wave Refraction using Fortran 63 
on the C.D.C. 1604 Computer 
Subroutine Title: FRICT^ 

Variables of Subroutine : 



HHQ.... Height of wave at point X, Y after bottom friction 

has been included. 

RLO Deep water wave length. 

,RLQ Wave length at point X, Y. 

< V CFF Friction factor between point n and n+1. 

^RCFF (1 - Kf) 

TRCFF Total of (1-Kj) terms to point n. 

* : TCFF; Total friction factor to point n. 



^Summary of Subroutine : 

The subroutine calculates the value of Kj and wave height reduced 
by bdttom friction (HHFR) for each point along wave ray. 



See Stouppe (1966) for discussion of remainder of program. 
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10 CALL KFUNCT(AA,FK;K) $ FKBAR=.5*(FK FKK) 00002520 

IF( 11-13)5,37,9 0000253 

37 FKKPP=FKBAR 000025A 

5 I F ( MAX-1 ) 7 , 7 , 9 0000255 

7 IF( IT-1 )20,20,9 0000256 
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IDENT SINHF SINHFOOl 

system SINHF002 

ENTRY SINHF SINHF003 

EXP DEC 709.089 SINHF004 

EXPl DEC -709.78 SINHF005 
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65 SINHF SINHF073 

end SINHF074 
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APPENDIX 2 



Sample Computer Input and Output 
Sample Input 



52 56 



A40. 


460. 


515. 


505. 


510. 


470, 


, 425. 


355, 


. 300. 


330. 


1 370. 


. 360. 


340. 


345. 


350, 


365. 


375. 


390, 


, 455. 


560, 


, 700, 


850. 


. 950, 


, 865. 


550. 


490. 


470. 


460. 


415. 


380, 


, 335. 


400, 


. 460. 


520. 


> 560, 


. 595. 


510. 


460. 


400, 


315. 


275. 


300, 


, 220. 


140, 


80. 


55. 


53, 


1 51. 


360. 


405. 


450. 


435. 


445. 


380, 


, 330. 


270, 


, 235. 


250, 


, 280. 


. 280. 


290. 


297. 


300. 


305. 


315. 


360, 


, 420. 


490, 


, 655. 


810. 


. 930. 


. 870. 


540. 


500. 


515, 


. 515. 


480. 


445, 


, 440. 


500, 


. 540. 


593, 


. 630, 


, 630. 


610. 


560. 


515, 


. 400. 


315 . 


275, 


, 234. 


190, 


, no. 


63, 


, 53, 


, 50, 


295. 


370. 


370, 


. 325. 


325. 


307, 


, 260. 


190, 


, 170. 


175, 


. 190, 


, 190. 


240. 


250. 


250, 


. 265. 


275. 


305, 


, 405. 


470, 


, 570. 


750, 


. 900, 


. 885, 


560. 


540. 


550, 


. 560. 


550. 


540, 


, 550. 


605, 


, 620. 


640, 


, 610, 


, 570. 


540. 


590. 


570, 


, 430. 


300. 


310, 


> 280. 


240, 


, 170. 


70, 


, 55, 


, 50. 


240. 


250. 


235. 


195. 


225. 


200. 


190. 


160, 


130, 


140. 


160. 


170. 


190. 


195. 


210. 


230. 


250. 


280. 


350. 


427, 


490. 


670, 


820, 


860. 


610. 


590. 


590. 


615. 


600. 


610. 


620. 


660 . 


650. 


620, 


560. 


505. 


480. 


520. 


560. 


515. 


405. 


370. 


320. 


260, 


190, 


80. 


60, 


49. 


150. 


195. 


200. 


110. 


125. 


120. 


115. 


95. 


80. 


98. 


105. 


no. 


162. 


170. 


180. 


205. 


210. 


210. 


290. 


380. 


455. 


580, 


740. 


785. 


660. 


650. 


640. 


660. 


660 . 


670. 


690. 


665, 


650, 


585, 


510. 


480. 


450. 


490. 


550. 


525. 


450. 


380. 


325. 


270. 


215. 


120. 


65, 


50. 


80. 


80* 


90. 


75. 


73. 


70. 


72. 


70. 


70. 


70. 


72, 


80, 


110. 


120. 


140. 


160. 


175. 


195. 


240. 


290. 


390. 


475. 


610. 


670, 


760. 


720. 


715. 


725. 


730. 


710. 


680. 


655. 


620. 


540. 


480. 


445 . 


440. 


410. 


550. 


490, 


415. 


370. 


325. 


285. 


235. 


170, 


85. 


65. 


61. 


60* 


60. 


58. 


60. 


60. 


61. 


63. 


60. 


60. 


65. 


65. 


80. 


85. 


100. 


120. 


130. 


120. 


160. 


2 30. 


280. 


380, 


460, 


540. 


770. 


750. 


735. 


715. 


665. 


635. 


620. 


610. 


570. 


485 . 


440, 


415. 


425. 


490. 


550. 


500. 


445. 


375. 


315. 


295, 


245. 


190. 


no. 


80. 


47. 


47. 


48. 


49. 


52. 


54. 


53. 


54. 


56, 


59, 


60. 


60. 


63. 


65. 


70. 


77. 


82. 


85. 


110. 


150. 


185. 


250. 


340, 


395, 


560. 


610. 


590. 


600. 


570. 


540. 


525. 


520. 


500. 


460. 


410. 


380. 


420. 


510. 


515. 


500. 


465. 


410. 


283. 


265. 


235. 


210. 


145. 


105'. 



4500. 


2 


, .01 


0, 


5.0 18,0 


10.0 


10 


30.0 


1500.0 


1500 


4500. 


2 


. .01 


6. 


5.0 18.0 


10.0 


10 


30.0 


1500.0 


1500 
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SAMPLE COMPUTER OUTPUT 
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SAMPLE COMPUTER OUTPUT 




• 001 OO'J 003 004 005 OOC 00? 006 



K-SCflLE - 1.00E+04 UNITS- IMCH. 
^ -SCALE ' 1.00E-fO4 UNI'^'^-'INCH. 



FAPRELL UflUE REFRACTION PROGRAM MONTER EY BAY 
"IDE LUL 0 ANGLE = 0 PERIOD = 18 



APPENDIX 3 



SUMMARY OF EQUATIONS 

Wave Speed (C) 

C |Z w. (^) . 

Ray Curvature (FK) 

FK = Z ^ . 

Coefficient of Refraction (K^ 

The value of ray separation {/S ) is calculated by solving the 
second order, non-linear differential equation: 



0 - + 

where 

p - - ^ ^ /) f ^ 

The above equation is solved by the finite difference method. This 
results in an equation for the Beta value at the nfl point in terms 
of the Beta values at the two previous points. The equation to be 
solved is then: 



where 



^ _ (pt> -Z)X3, 

' 2+pP 






D is the incremented distance 
p,q are as defined above 

'^1, ^2 Beta values of the two previous points. The 

The coefficient of refraction is calculated by the relation: 
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Direct Shoaling Coefficient (Xg) ; 



s.zsn 



I z. + Z S' ~ 

+ II. h ^ 



Wave Height Reduction Factor Due to Bottom Friction (Kf) : 



where 







A//. + I 



) 




-F H. 

^ L 
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